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Abstract 

T(3)-gauge model of defects based on the gauge Lagrangian quadratic in the 
gauge field strength is considered. The equilibrium equation of the medium is 
fulfilled by the double curl Kroner's ansatz for stresses. The problem of replication 
of the static edge dislocation along third axis is analysed under a special, though 
conventional, choice of this ansatz. The translational gauge equation is shown to 
constraint the functions parametrizing the ansatz (the stress functions) so that the 
resulting stress component 033 is not that of the edge defect. Another translational 
gauge equation with the double curl differential operator is shown to reproduce 
both the stress functions, as well as the stress tensors, of the standard edge and 
screw dislocations. Non-linear extension of the newly proposed translational gauge 
equation is given to correct the linear defect solutions in next orders. New gauge 
Lagrangian is suggested in the Hilbert-Einstein form. 
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1 INTRODUCTION 



Considerable attention has been paid in last years to various applications of the ISO (3) 
(IS 0(3) is the group of rigid body T(3)x> SO (3)) gauge model of defects in (continuous) 
solids proposed in [1] and, as a more elaborated version, [2]. For instance, the gauge 
theory of continuum damage in solids [3] , the gauge theory of elastic materials exhibiting a 
relaxation phenomenon [4] , the gauge theory of "plastically incompressible" elastic-plastic 
non-dissipative medium [5] have been considered, as well as the problem of electronic states 
in solids containing isolated defects [6, 7]. All these applications [3-5, 7] have been mainly 
concerned with the translational sector of the model [1, 2], i.e. with the Lagrangian L§ 
quadratic in the gauge field strength postulated in [1, 2] to govern the translational gauge 
field. 

The applications found look promising for the class of the T(3)-gauge models with 
the quadratic Lagrangians like C^. On the other hand, a conventional approach to the 
defects in solids, i.e. the defect theory [8] which is concerned with the Volterra-type 
solutions of incompatible elasticity theory, is widely acknowledged. The theory [8] deals 
with the singular dislocations and disclinations, as well as with their distributions, and it 
admits a number of reliable calculations in the field of mechanical problems of solid state 
physics [9] (for instance, and refs. therein). Besides, the use should be mentioned of the 
singular solutions [8] in the problem of sound propagation in elastic body with topological 
defects [10], and in such microscopic problem as particles diffusion in solid with randomly 
distributed defects [11]. For instance, the screw dislocation solution has been adapted 
both in [10] and [11]. 

Therefore it is seen that the conventional approach [8], as well as the non-conventional 
one [1, 2], are actively developed. However, with regard to the defect modelling in the 
framework of the Lagrangian T(3)-gauge approach, and namely with the quadratic 
the following problem lacks a sufficient attention by now: whether some of the solutions 
considered in [8] can be reproduced ('replicated', the terminology of [12]) in the formalism 
[1, 2]? The Refs.fl, 2] themselves give a confirmative answer about this both for the screw 
and edge dislocations (one should also refer to the original paper [12]). Moreover, the last 
statement has stimulated the attempt [7] to calculate the second order corrections to the 
screw dislocation in the formalism [1, 2]. 

Conventional theory of defects [8] can be naturally considered as the Abelian gauge 
model [13-15] with the additive gauge group iso(3) ~ M 6 [16] (iso(3) is the Lie algebra 
of IS 0(3)). In this case the motor [17] of the disclination and dislocation loop densities 
plays the role of the 6-component Abelian gauge potential, and the motor of the defect 
densities that of the Abelian gauge field strength [18]. It can presumably be concluded 
from [18] that when ISO (3) is gauged, a successful replication (as soon as the last can 
be questioned) of the defects [8] should be valid both for the translational and rotational 
ones. Unfortunately, any attempt to obtain the standard wedge and twist disclinations 
[8] in the formalism [1, 2] is seemingly absent. 

The present paper has originated as an attempt to understand, whether the Lagrangian 

[1, 2] (i.e., practically, the translational model [12]) leads to an extension of the classical 
picture [8] for (singular) dislocations, or it provides a model of its own significance with 
the solutions which should be interpreted on their own rights. 

Formally, the double curl Kroner's ansatz for stresses is used in [1, 2] to fulfil the 
equilibrium equation of the medium. In its turn, the translational gauge equation (the 
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gauge equation, for brevity) determines the stress functions which parametrize this ansatz 
since can be re-expressed by means of the constitutive relation 'strain-stress'. As to 
the screw defect, the gauge equation has been reduced in [1, 2] to the homogeneous 
Helmholtz equation instead of the standard Poisson equation for the stress function. In 
the axisymmetric case (8/8x3 = 0) the Helmholtz equation has provided the modified 
Bessel function K (p) as that solution, which coincides with the Prandtl stress function 
[19], i.e. behaves as log(l/p) + const, at p <C 1. Thus the screw defect is permissible in 
[1, 2]. 

In spite of the fact that the approach [12] to the screw dislocation is not contradictive, 
the situation with the edge one is more subtle. Generally speaking, the quadratic 
becomes inappropriate as soon as the stress functions of the edge dislocation (along third 
axis) are obtained in the same way as for the screw one. Practically, the T(3)-gauge equa- 
tion constraints the appropriate stress functions so that the resulting stress component 
o" 33 is not of the edge defect. 

An alternative form of the gauge equation is proposd in the present work, which 
eliminates the contradiction and, without any artificial effort, admits the linear solutions 
for the edge and screw dislocations. It is crucial that both the gauge equations, i.e. the 
newly proposed and that from [1, 2], being expressed through the general tensor of stress 
functions, pass each into other under the exchange v < — > is -1 , v is the Poisson ratio. This 
is just the point needed to avoid the problem found for 0^3, thought to keep the situation 
with the screw dislocation unaffected. 

The new gauge equation allows generalization, since can be written as an Einstein- 
type (non-linear) equation. This requires to modify the geometrical arena of the model [1, 
2], and suggests in the Hilbert-Einstein form. The stress tensor of non-linear elastic 
body plays the role of the source in the Einstein-type gauge equation. As a specialization 
of this source, the Murnaghan's constitutive relation of isotropic body is considered in 
the present paper. The new gauge equation should be appropriate to find out quadratic 
corrections to the linear solutions outside of the defect cores. 

The paper is written in five sections. Sec.l is introductory, Sec. 2 is concerned with the 
difficulties of the replication of the static edge dislocation, and a modified linear model 
(i.e. another translational gauge equation, in fact) is proposed in Sec. 3. Further, Sec. 4 
is devoted to the Lagrangian formulation of the corresponding non-linear model, and 
suggests in the Hilbert - Einstein form, while Sec. 5 concludes the paper. In what 
follows, it is assumed that reader can be referred to [1, 2] for motives of the ISO(3)- 
gauging, as well as for certain details and comments on the original formulation. The 
consideration is time independent, Greek indices run from to 3, Latin ones - from 1 
to 3, and repeated up and down indices mean summation. Besides, the Latin letters a, 
b, c, d,... denote curvilinear indices of a deformed configuration, and i, j, k, 1,... denote 
Cartesian frame of initial (undeformed) one. 
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2 T(3)-GAUGE EQUATION FROM THE QUAD- 
RATIC LAGRANGIAN AND RELATED DIFFI- 
CULTIES 

Following [1, 2] let us briefly remind the basic elements of the theory of the translational 
gauge field a a assuming that the gauge-rotational degrees of freedom of the whole ISO(3)- 
model are "frozen". The translational gauge field <p a a plays "compensating" role as the 
Abelian group T(3) acts by non-homogeneous shifts on the current configuration variable 
£ a , and the corresponding compensated derivative has the form B a a = d a £ a +(f) a a , where d a 
are the partial derivatives d/dx a , and xo implies time. The current configuration £, a (x, xo) 
can be written in the form £ a = 5f x % + u a (x,x ) with respect to a Cartesian coordinate 
system so that x implies an initial (undeformed) configuration, and u is displacement field. 
It should be noticed that when Abelian non-compact group like T(3) ~ M 3 is gauged, 
appearance of B a a instead of the purely gradient quantity d a ^ a implies that the latter 
ceases to be adequate variable, because non-integrable contributions become essential. 

As far as the problem of replication of the static edge dislocation is in focus of the 
present paper (see [8] for the standard results), it will be assumed that do = 0, (f) a = 
so that the governing equations are: 

<W = o, (i.i) 

d 3 - d^j) = dfiK - d^drfj) = (2 a )-V- (1-2) 

The source field a a l (1.2) is given by 

a; = B aj (\5^E k k + 2pE*) , (2) 

A and /j, are the Lame constants, and the strain field E^ is defined as 2Eij = B k i B^j — 5^ 
(see Sec. 3 on the usage of up and down indices). Notice that the stress field cx y is just 
given by the brackets in the R.H.S. of (2). It is seen that the equilibrium Eq. (1.1) ensures 
integrability of (1.2). The parameter s in the R.H.S. of (1.2) (the "coupling" constant) 
is due to the following choice of the spatial part of the quadratic gauge-translational 
Lagrangian: 

£ = (-23)0^^0^, (3) 

where the square brackets imply antisymmetrization. It has been stated in [1, 2] that 
the translational gauge equations (1.2) "are what replace the compatibility conditions of 
linear elasticity theory" . 

The Refs.fl, 2] propose the following approach to the dislocation problem. The purely 
integrable contribution to the state due to Uf, is zero. Therefore, Eq.(2) can be linearized 
in <f> ib , since B ib is 5 ib + (f) ib , as follows: 

Oij = 2/i0 ( y) + X5 i:j (f) kk , (4) 

where the brackets mean symmetrization (f>^ = (1/2) (0^- + (pji) (in linear case we do 
not distinguish up and down indices). Further, we impose the conditions 0i2 = 4>2i, 
013 = 03i = 0, 023 = 032 = 0, and suppose 83 = to adjust the axial orientation. In this 
case only a part of Eqs.(1.2) survives: 

a| 2 0n - <9r20i2 = (2s) 

#1! 022 ~ <9l2012 = (2s)" 1 0- 22 , 
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A033 = (2s)-V 33 , (6) 

where A = <9 2 x + d 22 . It can be realized that the choice of <pij in the symmetric form 
<t>ij — 0(ij) should provide us the appropriate solution of (1.1), (5), (6) as a small elastic 
strain tensor, while Eq.(4) gets the status of the constitutive relation of isotropic body 
(the Hooke law). 

The way to handle the system (1.1), (5), (6) proposed in [12] is as follows: to postulate 
the ansatz 

/ d 2 22 f -dl 2 f o\ 

VT X ° = -dnf %J , (7) 
V Op) 

which fulfils (1.1), while Eqs.(5), (6) enable to specify the unknown functions / = f(xi,x 2 ) 
and p = p(xi,x 2 ) provided the use of the linear law (4) is made. 

As far as we are interested in the edge dislocation along the third axis Ox 3 , the 
component 033 (i.e. the strain (33)-component) is expected to become zero, at least as 
a limit of an "extended" solution. As 033 = u(an + <J 22 ), where v = A/2(A + p) is the 
Poisson ratio, implies vanishing of the strain (33)-component, the constraint uAf = p has 
to be expected in that limit. The relation (7) with p = uAf but, formally, with opposite 
sign at / is nothing but completely the standard ansatz of the theory of dislocations, 
where / is called Airy's stress function [19]. This function / has been found in [19] as the 
biharmonic potential <92(p 2 logp) (or <9i(p 2 logp)), p 2 — x\ + x 2 ,, and it enables to obtain 
via (7) all the stress tensor components of the edge defect. 

Let us re-consider the solution [1, 2] of the problem of the edge dislocation. Clearly, 
Eq.(l.l) is respected by a (7), whereas Eqs.(5) result in 

(l-a)A/ -ap= K 2 f, (8) 

and (6) leads to 

(1 — a)Ap — a A Af = n 2 p , (9) 

where a = A(3A + 2p) _1 and k 2 = p/s. Let us exclude p from (9) with the help of (8). 
We assume here the limit k — > proposed in [12] to restore the classical situation [8]. In 
this case AAf = arises to define the limiting form of /, while (8) provides 

P = — A/=-A/ (10) 
a v 

as the limiting form of p. 

The point is that the fourth order differential equation, which appears to govern / at 
arbitrary k, looks troublesome for analytical solution. It is why, in view of AAf = 0, the 
authors have been forced to "guess" that / can be found so that the needed Airy's function 
behaviour appears for it at k — > 0. But since AAf = is to indicate at the correct /, 
just the same reasons (i.e. linearity of the substitution of p) imply that Eq.(10) defines 
the limiting form of p. However, the Eq.(10) contradicts the constraint demonstrated 
above. Thus, in spite of the satisfactory equation for /, the correct stress component a 33 
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of the edge dislocation would not appear neither at finite k nor at k — > 0. Notice that in 
the Ref.[l] it has been recognized only for finite k that the new <t 33 is different from the 
standard one. 

Otherwise, let us try "rigidity" of (8), (9) against the constraint p = vAf by substi- 
tuting it to them: 

(l-i/) A/ = n 2 f, 
= n 2 Af . 

These equations can hardly be fulfilled with nontrivial / at finite k, and, moreover, 
nontrivial solutions of AAf = with the analytical behaviour we are interested in would 
not appear even at k — > 0. 

Therefore, the joint use of Eqs.(5), (6) and (7) to obtain the edge dislocation's stress 
function leads to the following conclusion. The resulting Eqs.(8) and (9) constraint the 
stress functions / and p so that it is impossible to get / as the biharmonic potential, and 
to fulfil at the same time p = vAf. One gets either only a 33 is incorrect due to (10), or / 
is not the Airy function at all. Besides, there is a little hope that an improved / can be 
found at finite k so that the additional contribution (— 1/a) lim^o (^ 2 /) from (8) to the 
R.H.S. of (10) could restore the desired p = vAf. 

With the purpose of illustration, let us obtain the limiting form of the matrix </> by 
means of (4) inverted and (7). We shall substitute the known / = (— C/2)(9 2 (p 2 logp), 
where C = (b/2n)(l — f) -1 , but p will be evaluated accordingly to (10). Then we get: 



4> C 



2p 2 



/ x 2 (l-2f) Xl {l-2f) \ 

x 1 {l-2f) x 2 (l + 2^|) 
V ' 2(l-±)x 2 ) 



(11) 



Integrating the equations d±u = 0n, d 2 v = (p 22 , d 2 u + d\V = 20 12 with respect to the new 
fields u, v, we obtain: 

v + iu = v d + iu d + -^-\og(x 1 - ix 2 ), (12) 

where i = \J — 1 , and u d , v d are the edge dislocation's displacements in the plane perpen- 
dicular to Oxs (the Burgers vector is in xi-direction) [8]. Equation (12) is to express that 
namely the logarithmic term responsible for the "closure failure" just characterizing the 
defect does not enter to u, v, because v d + iu d contains (—b/2n) log(xi — ix 2 ). Besides, 
033 (11) looks embarrassantly to consider it as the corresponding strain component of the 
plane problem in question. In spite of the statement by [1, 2], it is plausible to conclude 
that Eq.(ll), as the limiting form of possible solution to (5), (6), does not imply at all 
replication of the edge dislocation at k — > 0. 

It can occur that additional care with k — * is needed to make rigorous statement 
about absence (or presence) of the edge dislocation solution for the Eqs.(5), (6). The 
formal implications of the use of Eqs.(4), (5)-(7) have only been pointed out here. In 
the next section a modification of the situation, i.e. a way to keep the ansatz (7), to 
exchange the parameters a and 1 — a, and thus to avoid the unpleasant p = z/~ 1 A/, will 
be demonstrated which admits the standard dislocations in the formal limit k — > 0. 
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3 THE DOUBLE CURL T(3)-GAUGE EQUATION 



The equilibrium equation (1.1) can be satisfied identically provided the stress field a is 
chosen as double curl of a twice differentiable symmetric tensor potential x which is called 
[19] the tensor field of second order stress functions: 

a i3 = (incxh = -6 ikl ei™d 2 km xin. (13) 

Both the particular ansatz proposed in [1, 2] to obtain the edge and screw dislocations 
appear as specializations of the general representation (13) (though (13) itself is not 
stressed in [1, 2]). 

The so-called stress function method based on (13) has been developed in [20, 21] to 
approach to non-linear dislocation problems. For the internal stress problem (isotropic 
case) the representation (13) has also been discussed by Kroner in [19]. The concept of a 
torsion-free stress space and relationship to it of the Eq.(13) have been considered in [19]. 
In [22] a stress-strain duality has been presented which relates the double curl ansatz (13) 
to a similiarity in Riemannian descriptions of kinematics and statics of mechanical state 
of solid. It has been exploited in [22] that the tensor of secon order stress functions and 
the stress tensor could play the role of the metric and the Einstein tensors, accordingly, 
of a non-Euclidean stress space. 

Therefore Eq.(7) looks more fundamentally than Eqs.(5), (6) with regard to the un- 
satisfactory implications of p (10). On the other hand, it is just the choice of the master 
equation (1.2) for the translational field which would require a justification. It turns out 
that the most direct way to get the function p consistent with the requirement 33 = 0, 
i.e. (T33 = v{(Jxy + 022), at k — > 0, is to choose the L.H.S. of Eq.(1.2) (linear consideration) 
as follows: 

(inc = - e ikl Od») = (2s)-^« > ( 14 ) 

where superscript Sym implies tensor symmetrized, e lkl is the permutation symbol, and o^- 
is given by (4). When the R.H.S. of (14) is zero, it looks like a compatibility equation of 
linearized elasticity theory provided 0^ is considered as distortion. Equation (14) replaces 
the compatibility condition in the sense that its R.H.S. is non-trivial while the L.H.S. looks 
conventionally. In the next section Eq.(14) will be viewed as the Einstein-type one. 

Let us accept (14) instead of the linearized (1.2) in the T(3)-gauge model of defects. 
However, it is more appropriate to rewrite (14) through the unknown second order stress 
functions %. To this end one should express (p Sym using both the Hooke law (4) inverted 
and (13), and, after this, substitute it into the L.H.S. of (14). The stress field cr (13) takes 
its place in the R.H.S. of (14). Finally, the fourth order differential equation appears: 

A< s >A< s > Xy + aD y A< 3 > X + {{l-a)d* + a - A< s > (fe + d f kXik ) = 

= k 2 (A (3) Xij + DijX + SijdkiXki - d.lxjk - df k Xik) , (15) 

where x implies trace of a is defined by Eq.(10), k is defined in (8), (9), and we have 
denoted the differential operators A^ 3 ^ = 5ij <9J, and D^- = — 5ij A^ 3 \ In its turn, the 
Eq.(15) can be considerably simplified if we replace x by another symmetric potential x' 
as follows: 
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where x' fulfils = 0. Thus we get: 



To clarify the situation let us repeat the procedure leading to (15) for the Eq.(1.2) as 
well, though assuming that the latter is written in terms of <p Sym . It is curious that thus 
resulting equation is completely similiar to (15) except only one thing: a and 1 — a are 
just exchanged. In other terms, v and are exchanged. Therefore it can be guessed 
that the problem related to p (10) should disappear when (14) is used instead of (1.2). 

It has to be noticed that linearizing Eq.(1.2) one gets another, apart from the problem 
of p, unpleasant thing: the source cr^ (4) in its R.H.S. is symmetric in the indices, while 
the L.H.S. is obviously not. It is why in [1, 2] <pij has been additionally assumed symmetric 
to perform the concrete calculations. On the contrary, the L.H.S. in (14) is symmetric, 
and the antisymmetric part of the translational gauge field automatically drops out, but 
must be taken into account to restore, say, distortion from strain field. 

In order to specialize (15) to the edge dislocation case, we introduce two functions p 
and / as follows: 

/ = X33, P = -d£ 2 Xll ~ 9^X22 + 2012X12, 

while other Xij are zero. The plane problem is adjusted by <9 3 = 0, and thus A^ 3 ^ becomes 
A (see Sec. 2 for A). It can be verified that now six equations (15) are reduced to 

aAf + (l-a)p = K 2 f, (16) 

(1 - a) AAf + aAp = k 2 p . (17) 

More precisely, Eq.(16) appears instead of the three equations d 2 ± A = 0, d 2 2 A = 0, and 
d 22 A = where A = (1 — a)p + aAf — k 2 f . It is seen that changing / to — / and 

a < > 1 — a one obtains (8), (9) from (16), (17), accordingly. 

Therefore Eq.(14) has provided us with the remarkable opportunity to convert the 
embarrassing ratio (1 — a) /a in (10) into a /{I — a). Finally, the resulting function p 
(16) gets around the obstacle discussed in the previous section. Notice that choosing 
/ = — <9iX23 + 02X31 one can deduce from (15) a single equation 

Af = n 2 f, (18) 

which is the same as in [1, 2] for the screw dislocation ansatz. It is clear that this 
coincidence is because (18) does not contain the Poisson ratio v — a/(l — a). Solutions 
to (16)— (18) (and, generally, to (15) ) should be called "modified stress functions" to 
distinguish them from the classical harmonic and biharmonic potentials. 
As the Eq.(16) defines p, the second equation governing / gets the form: 

(a-m 2 ) (a + a/- 2 )/ = o, ^ 2 = 7Y^- (19) 

For correspondence with [1] we use M 2 instead of our k 2 when considering (19) and its 
solution. Equation (19) simply differs from the corresponding Eq. (4.6.27) in [1]: only 
2a — 1 and 1 — 2a are interchanged under a < — > 1 — a. 
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As 82 commutes with A, one can be concerned with / in the form c^/^p), p = \x\, so 
that h also respects (19) with A reduced to 

-4- (4-) 

p dp \ dp ) 

Equation (19) has the Bessel and Neumann functions Jo(J\fp) and Y (Afp), and the mod- 
ified Bessel functions Io(Aip), Ko(Aip) as four angle-independent basic solutions [23]. 
These basic solutions can be combined so that their combination is decreased at infinity. 
Eventually, the following solution can be obtained: 

(20) 

Hp) = log^ Jo(A/» - \y q (M P ) - K (Mp). 

At Aip,Afp <C 1, the solution / (20) results in the standard biharmonic potential, since 
T ~ M 2 /(2(l — v))p 2 log p. Notice, that the normalization of T(p) is chosen so that 
the stresses are a n = — <9f 2 /, etc., but not, accordingly to (7), a n = —pd 2 2 f : etc.. The 
limit A4p,J\fp <C 1 implies that p is finite but not zero while p/s — > 0. In the opposite 
case A4p,J\fp ^> 1, the solution (20) becomes a linear combination of X2 p~ 3 ^ 2 sin(jVp) and 
%2 p~ 3 ^ 2 cos(jVp), i.e. it is O ((A/"p) _1 ^ 2 ) • Therefore, we are expecting all the entries of the 
stress matrix cr at A4p,Mp <C 1 to be those of the edge dislocation, because p = uA(—f) 
by k — > in (16). Recall that the solution to (18) K , which enables the screw dislocation, 
possesses fast and monotonic decreasing as p~ 1//2 exp(— up). In other words, this solution 
is really 'short-ranged' in comparison with the corresponding classical one. 

Let us also mention the Ref.[39], where a "mass" of defects is discussed as an im- 
plication of the translational model [1] (specifically, as implication of the corresponding 
time-dependent T(3)-gauge equation). This is because the Klein-Gordon equation is pos- 
sible with the mass k 2 for the translational gauge field in the Lorentz-like gauge. The 
present section demonstrates, that such "mass" effect should be rather traced to Eq.(15) 
for the modified stress functions. The point is that namely (14), but not (1.2), seems to 
be related to the defects (at least, to the statics of the conventional ones). Moreover, it is 
just the Eq.(15), which implies the conventional, unboundly increasing, stress functions 
to become the limited modified ones just due to the "mass terms" k 2 in (18) and (19). 

However, in Sec. 4. 3, and 5, it will be discussed, that only the replication demonstrated 
above (i.e. the regime up <^ 1) is of main interest in the linearized version of the Einstein- 
type consideration. Thus, the L.H.S. of (14) helps to avoid the problem with p (10). Both 
the sides of (14) admit non-linear generalizations, and therefore it is suggestive to put 
(14) into a more general form which can be derived in Lagrangian approach. 

4 THE EINSTEIN - TYPE EQUATION AND ITS 
LAGRANGIAN DERIVATION 
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4.1 The geometric preliminaries 

Before to proceed with the Lagrangian derivation of Eq.(14), let us briefly present the 
geometrical apparatus which will underly our main construction. To make the discussion 
reasonably compact, it will be assumed that one can be referred, say, to [24] for basics 
of geometry of the Riemann and Riemann-Cartan spaces in the form accomodated to 
describe defects. Other useful references can also be found in [25]. As to gauging of the 
important for us group ISO(3), a formally close gauging of the Poincare group (which also 
is a semi-direct product of translation and pseudo-orthogonal rotation groups though of 
4-dimensional Minkowskian space-time) have already been extensively developed in the 
realm of gravitational physics [26]. One should consult with [26, 27] for a similiar, though 
much more elaborated machinery. 

Although we are restricted to the T(3)-gauging, it is more appropriate to admit, just 
for a moment, a more general framework of the 7S'0(3)-gauging. Here the couple of the 
Cartan structure equations 

(21) 

T% = d t B a 3 - d 3 B% + co^B^ - u/^-BV 

appears as one of the basic differential-geometric relations. Four tensors which enter the 
Eqs.(21) are components of the differential forms which define the geometry of the so- 
called Riemann-Cartan spaces: H a bi j an d T a %ji which are antisymmetric in determine 
the curvature and torsion 2-forms, accordingly H a b and T a . The coefficients uo\^ and B a i 
define the connection 1-form u a b and the co-frame 1-form B a , respectively. The couple 
of indices a, b demonstrate TZ a b , oj a b as elements of the Lie algebra of the group 5*0(3), 
while T a , B a belong to the Lie algebra of the group T(3). All the fields are considered in 
open domain of a 3-dimensional manifold with Euclidean signature. 

It is known that gauging ISO (3) one gets the iso (3)-valued connection and curvature 
differential forms, both of which are split into "linear" and "translational" parts owing to 
the semi-direct sum structure of iso (3) [27, 28]. However, the "translational" curvature is 
transformed non-covariantly under the IS 0(3) gauge transformation. In order to find out 
the quantity which is transformed (gauge-)covariantly, it is proposed to use the auxiliary 
field £ a [27-29]. As the result, the Cartan structure equations in the ISO (3) gauging 
acquire the form (21) provided the coefficients B a i have the following structure: 

B a i = r i + d i e + u\ i e. (22) 

In (22) 0J a h ^ and <p a i are just the "linear" and "translational" parts of the iso(3)-valued 
connection 1-form. In other terminology, u a b i and a j are the gauge potentials which are 
due to non- homogeneous action of IS0(3). Going further, apart from it is also ap- 
propriate to define their reciprocals B b J as follows: B a i B b t = 8 b , B a i Bj = 5j. The point 
is that the components B b J define an orthonormal triad Bb = B b dj (more exactly, B b J are 
transitions between the coordinate basis in the tangent space {dj} and the orthonormal 
basis {B b }), whereas B a i provide the dual basis of 1-forms B a = B^dx 1 . 

It is crucial that the auxiliary field £ a (22) has been identified in [1, 2] as the deformed 
configuration variable (see [28] for other problem-motivated interpretations of that special 
field) as follows: 

£ : x i _ C(r') , r 1 : C — A?) . (23) 
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This equally means that the group indices a, b, c in (21) get the "material" meaning since 
the sets {£ a } label points in deformed configuration. The coordinate indices i,j, k in (21) 
get the status of the Cartesian ones in initial (undeformed) configuration. In particular, 
both B a i and £> a l are useful to pass from one set of indices to another. 

In [18] the ISO (3) Cartan's equations (21) have been reduced to the conventional 
expression of the motor of dislocation and disclination densities through the motor of 
the defect loop densities. It is essential that the field £ a which is x a + u a has enabled 
this truncation so that sequence of Schaefer's exterior differentiations inherent to [8] has 
naturally appeared. Besides, interpretation of the /S'0(3)-connection in terms of the 
defect loop densities has become possible. 

Regarding the material interpretation of {x 1 } and {£ a }> let us recall how the Green 
deformation tensor and the corresponding Lagrangian strain tensor Eij appear in 
the conventional elasticity theory [30-32]. Indeed, the length element in a deformed 
configuration can be written in the Euclidean form ds 2 = i] ab d^ a dC, b (i] ab is flat metric) as 
well. At the same time, the configuration {£ a } can be considered in terms of the initial one 
by means of the maps (23). In this case the length element becomes expressed through 
the initial coordinates:^ 2 = r] ab di^ a dj^ b dx l dx j . Thus we obtain the corresponding non- 
Euclidean metric (the Green deformation tensor) g^ = i] ab d^dj^, while the Lagrangian 
strain tensor is 

1Eij = gij - Vab5?Sj = 9ij - Vij ■ (24) 
Following the close analogy, we define the Green deformation tensor as 

9ij = rj ab B\ B\ , (25) 

with B a i (22), while the strain tensor is given by (24). The metric and its inverse can 
be used for raising and lowering the indices while r\ ab - to handle analogously a, b, c. 

Now it is time to stress that the present work is concerned with the Eq.(14) to avoid 
the problem discussed in the Sec. 2. The given section is to extend straightforwardly 
(14) as an Einstein-type equation in the context of torsion-free Riemannian geometry. 
Therefore, given the metric g^ is (25), a unique linear connection without torsion can 
be associated with it so that g^ is covariantly constant. The corresponding condition 
^iQjk = implies the conventional expression for such T^-: 

21* = g kl (d i9jl + d m - d l9ij ) . 

The connection can be related with the S'0(3)-connection uj\^ by means of the relation 

^\ t = B b l (^ l B\-d t B a l ), 

which respects the covariant Constance of g^, and, further, allows to obtain from H a bti j 
(22) the Riemann-Christoffel curvature tensor R k mij : 

b,ij — Bb k mij > 

T>k a,-pk f) V^ -I- F n F^ F n 

mij u i jm j im ' in jm jn im " 

As the last step, we define the scalar curvature R, which is an important geometric 
invariant as follows: R = R/, R^ = R k ikj - The scalar R is just what we need to derive 
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(14) in the Lagrangian approach. From now on, we "switch off" the gauge-rotational 
degree of freedom uo ah ^ in B a i . 

The geometric presentation above is inevitably sketchy, since many appropriate things, 
like definition of principle bundle of affine (linear) frames, definition of associated vector 
bundle, gauge transformation rules, etc., are omitted. Necessary matter can be picked 
up from the literature cited, though more detailed and accurate geometric background of 
the non-Abelian model would require a separate presentation. However, as to the IS 0(3) 
gauge model [1, 2], the idea to pass from the "triad" representation (21) to the purely 
"coordinate" one, in order to discuss replication of the conventional dislocations, belongs, 
seemingly, to this paper. 

4.2 The Lagrangian derivation 

It is well known that the six compatibility equations of elasticity are, in fact, vanishing 
conditions of the six independent components of the 3-dimensional Riemann - Christoffcl 
tensor R l jki- Due to 3-dimensionality, one can equally use the second rank Einstein tensor 
G ij instead of R*^ as follows: G ij = (1/ A)e ikl e jmn R k i mn , where e ijk = y/~ge ijk and g is 
det(gij). Linearizing the Riemann-Christoffel tensor one gets the double curl expression 
like the L.H.S. of (14). 

Therefore, the Eq.(14) can be considered as a weak field approximation of an Einstein- 
type equation. Indeed, the L.H.S. of (14) looks like linearization of G' lj , while the stress 
tensor in its R.H.S. can acquire higher powers of E^, as a non-linear constitutive relation 
'stress-strain'. In other words, the situation is reminiscent to a gravitational equation 
which relates the Einstein tensor (though 4-dimensional) to a matter energy-momentum 
tensor as the source. Therefore both the sides of (14) can be extended, and the resulting 
equation admits a Lagrangian derivation. 

Now we can proceed with the Lagrangian approach. First of all, we postulate the 
Hilbert-Einstein Lagrangian density which is responsible for the translational field: 

— Cs = s R , 
b 

where b = <iet(-B a j), and R is the scalar curvature. Variation of takes the form (up to 
terms irrelevant by the Stokes theorem): 

I M± = 2s ( Rj - — Bj) = 2s G* , (26) 
b 5B% V a 2 a ) a ' K ' 

where G a * is the Einstein tensor. 

Let us obtain a in the R.H.S. of the Einstein-type equation. The field £ a is an im- 
portant constituent of more general /5'0(3)-formalism, and it should be governed by an 
appropriate Lagrangian C^. Our consideration is static, and therefore we shall choose 
(—l/b)C^ in the form of potential energy W of isotropic non-linear elastic continuum 
(practically, in the so-called Murnaghan's form). Given the field B a i has the T(3)-invariant 
form d£ a + a j, variation of the Lagrangian gets the form: 

— — = - £ J = - aj - B n l W . (27) 

b 5B a i a a a V J 
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In (27), by definition, the stress field a lj is given by SW/SEij 
in S a * is due to variation of b. 

Putting together (26) and (27), we obtain the equation 



G; = (2s) -1 E 



a i 



a l K and the second term 



(28) 



which generalizes (14) in the sense explained above. Rejecting in (28) higher terms and 
assuming coincidence, in the leading order, of its source with (4) (weak field approxima- 
tion), we obtain (14). Tending s to infinity, we just recover the general compatibility 
equation G^ = 0. As far as the Einstein tensor is "covariantly conserved" by the appro- 
priate Bianchi identity [24], the equation VjS a * = appears to govern the source tensor. 
For the first time in the context of gauge dislocations, an equation similiar to (28) has 
been obtained by variational approach in [33], but without any further elaboration. 

For definitness, let us specialize the potential energy W taking it in the form proposed 
by Murnaghan [34]: 

W = ^ll{E) + fih(IE 2 ) + ^ll{E) + v 2 h{E)h{E 2 ) + ^u 3 h(lE 3 ), (29) 
2 6 3 

where IE implies the gauge invariant strain tensor Eij (24), invariant function ii(...) 
implies trace of the appropriate tensor argument, while A and \i are the Lame constants 
of second, and v\ 2 3 - of third orders. The stress tensor acquires the form 



\h(E) + -±It(E) + v 2 h{E 2 ) +2E 13 [fi + v 2 h{E)\ + 4u 3 E 



pik 



El, 



(30) 



E k k and h(E 2 ) = E ki E kl . When u lj2i3 are zero, we obtain (2), while a. 



■ij 



where I±(E) 
coincides with (4) in the weak field limit. 

The problem of corrections to the linear dislocation solutions has been actively inves- 
tigated by various methods. For instance, the Refs.[20, 21, 35, 36], as well as a review of 
them in [31], are to be mentioned here. The Refs. [31, 32, 37] are also useful as the sources 
of information about others approaches to the non-linear dislocation problems. Moreover, 
in [38] second order corrections have been obtained for a wedge disclination solution. All 
the calculations in the refs. mentioned have been done assuming the constitutive relation 
in the Murnaghan's form to account for all the quadratic contributions properly. The 
Ref.[7] is concerned with the non- conventional approach to the screw dislocation [1, 2]. 
However the constants ^1,2,3 are taken zero in [7]. 



4.3 Vanishing torsion ? 

Before to conclude the paper, the Eq.(14) has been proposed in the Sec. 3 to avoid the 
contradiction between p (10) and 033 = v(o\\ + 022) for the edge dislocation. In their 
turn, Sec. 4. 2, 4.3, are to extend (14) to the Eq.(28), which is concerned with the Riemann- 
Christoffel picture. Therefore, it is appropriate to check, whether the solutions considered 
above are consistent with the requirement of zero torsion. Obviously, the torsion can not 
be asked to vanish everywhere, because thus defects will drop out. It seems sufficient, 
since our description should approximate a more adequate Riemann-Cartan situation, to 

ask about the torsion which is Z6rO, clt l6cLSt, dbt Kp <§C Rcore P -^exterior OO. 
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First, it is appropriate to remind the situation [12] with the screw defect, which is 
more transparent. It is known how to calculate the torsion tensor [24] (see also (21)), and 
so we obtain from (4) and (7) for its single non-zero component: 

T\ l2 = j [2(8^2 + c^) - A/] , (31) 

where the modified stress function / fulfills (18). There are two auxiliary functions iQ\ 
and 0J2 in (31) to account for the antisymmetric part of distortion: 

/ d 2 f-2co 1 

= - -drf + 2002 

1 \d 2 f + 2oor -dj - 2oo 2 

In the (anti-)plane problem we choose cu 1 , cj 2 so that the distortion components 3i and 
032 are zero (for the defect along Ox 3 ) [8]. Then (31) is reduced to 

T 3 , 12 = -A/. (32) 

In the axisymmetric case, the solution to (18) can be written as / = (b/2ir)K (Kp) [1, 
2]. At np <C 1, this solution results approximately in (b/2n) log(l/p) + const, i.e. in the 
Prandtl's stress function, while the corresponding value of T 3 5 12 (32) is zero, since we get 
bS(xi)S(x2) at p 7^ 0. At up ^> 1, the torsion T 3 5 i 2 approximately vanishes by (18), (32), 
since the Bessel function decays exponentially. Provided the entire / is used to consider 
the torsion component at up <C 1, one gets: 

T 3 5 i2 ~ - -^-k 2 (const - log(/tp)) . 

The latter still vanishes because p is limited R core < P < Rexterior an d k — > 0. However, 
T 3 12 (32) is not zero in the intermediate region up ~ 1, thus violating our Riemannian 
interpretation. 

Is it possible to make T 3 12 zero by an appropriate choice of u>i, oj 2 ? Formally, it can 
be done relaxing the requirement 03i = 03 2 = 0. But therefore the framework of the 
(anti-)plane problem is left in favour of a 3-dimensional consideration. Notice, that the 
torsion components have not been considered in [1, 2] for the screw dislocation as we just 
did. This has influenced, for instance, the Burgers vector component b 3 found in [1] as 
b/2 instead of b owing to the neglection of 3 i = 32 = at k — > 0. 

Now let us turn to the edge dislocation which is more complicated. Here the torsion 
components are: 



T\ 12 = ]- \2d 1 oj + (l-u)d 2 Af + un 2 d 2 f 



r 2 



12 



2 

1 
~2 



(33) 

2d 2 u - (1 - ^)diA/ - vK 2 d x f \ , 



where u is to account for the antisymmetric part in the distortion component 0i 2 (the 
normalization of / corresponds to (7)). The dependence on k is present in (33), because 
p (16) has been used. When /, T 1 12 , and T 2 12 are fixed, the function uj can be found 
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from (33) provided, say, the integrability condition (8182 — 8281)00 = is implemented, 
i.e. the equation 

(l-f)AA/ + uK 2 Af = 2(d 2 T\ l2 - d.T 2 ^) , (34) 

holds. 

In the formal limit k — > the choice 

b X\ 

ensures that (33) can be fulfilled with 

T\ 12 = b5(x 1 )5(x 2 ), T 2 , 12 = 0, (35) 

and / = (6/47r)(l — z/) _1 <9 2 (p 2 logp). The torsion, Eq.(35), is given, the integrability 
Eq.(34) acquires the form at k — > 0: 

(l-i/)AA/ = 2bS(x 1 )d 2 S(x2), 

thus prescribing the limiting form of the modified stress function. 

Taking into account the explicit solution / (20), and using (33), (34), we can also 
calculate the torsion components T 1 i i2, T 2 5 i2 at arbitrary k. Indeed, let us take T 2 ,12 = 
to express 00. Then we obtain T 1 i i 2 : 

T 1 12 = -7- f(l - v) AAF + vk 2 AT) = 

47T V ' 

where the Eq.(19) has been used to exclude AA. Using T{p) (20) explicitly, we obtain 
further: 

T\ 12 = ±J^-((1 + V )jr {p) + 2uK (Mp)) • (36) 

47T 1 — V 

At Mp,Afp 1, the R.H.S. of (36) is also vanishing due to the same reasons, as it 
happens for T 3 t i2 (32) in the analogous situation. At Aip,Afp ^> 1, T 1 i i2 decays as 
linear combination of p -1 / 2 sin(jVp) and p -1 / 2 cos(A/"p), i.e. not so fast as T 3 12 (32) does. 

Is it possible to make 12 , T 2 12 zero by adjusting ul The integrability condition 
for uo can be written, without the use of (16), as follows: 

(1 - a) A Af + aAp = 

Then, formally, at finite s the parameter p, i.e. the stress component 033, is zero by (7), 
(17). Therefore the strain component -E33 is 

£33 = + °22) = |A/. 

If E 33 is still zero (as at k — > 0), the function / is also zero by (16). Otherwise the 
two-dimensional consideration is no longer sufficient. 

Let us sum up this section. Here we have obtained the torsion components expressed 
through the stress functions. The results of the Sec. 3 have been used to establish that 
there exists a single limiting situation, i.e. the limit R core < p < R ex terior, « — ► 0, when 
the modified stress functions can be replaced approximately by the conventional ones, 
and the torsion becomes zero. This would imply that the Einsteinian interpretation of 
the Eq.(14) is valid only in that limit, while the modified stress functions themselves for 
others k, p should be discarded at the present stage. 
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5 DISCUSSION 



Thus, in the present paper we have started with the demonstration of the fact that the 
conventional linear solution for the edge dislocation does not fit in the gauge model [1, 
2] (specifically, its translational sector). The Kroner ansatz for stresses is used in [1, 
2] to fulfil identically the equilibrium equations of the medium. It is essential that the 
corresponding parametrizing functions, i.e. the (modified) stress functions, are to be found 
from the translational gauge equation. In its turn, the form of the translational gauge 
equation [1, 2] is dictated by the Lagrangian quadratic in the T(3)-gauge field strength. 
As a special example, in [1, 2] the Kroner ansatz has been chosen in that form [19], 
which leads in conventional theory to the Airy's stress function of the edge dislocation 
along third axis. The present paper displays that in this case the gauge equation [1, 2] 
constraints the modified stress functions so that the resulting stress component 033 arises 
as z/~ 1 (cr 11 + (T22) instead of ^(o"n + 022) at k — > 0, i.e. o" 33 is not of the edge dislocation. 

The ansatz itself looks more valuably for the problem in question, and therefore the 
contradiction found can be eliminated by another choice of the master translational gauge 
equation. Namely, the use of the double curl differential operator in its L.H.S. leads to 
the replacement a < — ► 1 — a, or, equivalently, v < — > v^ 1 . More precisely, both the 
linear gauge equations in question (our double curl Eq.(14)and that from [1, 2]), provided 
they are written through the general tensor of stress functions, pass each into other just 
simply by inverting the Poisson ratio v. Eventually, both the standard linear solutions for 
dislocations are possible with the new gauge equation, though in the formal limit Kp <C 1 
(k — > 0, p is finite). 

From a more formal point of view on Lagrangian field-theoretical modelling of con- 
ventional (and non-conventional) defects, the newly proposed Eq.(14) looks suitable for 
the replication of the dislocation solutions. Therefore, it is worth to be realized that the 
double curl in the L.H.S. of (14) can be considered not only formally, as a nice trick, 
but also as linearization of a three-dimensional Einstein tensor. Moreover, a non-linear 
constitutive relation 'stress-strain' can be thought of instead of the R.H.S. of (14). Thus, 
(14) can be treated as a weak field limit of a non-linear Einstein- type Eq.(29), which is the 
most straightforward extension of (14). The Eq.(29) uses the framework of torsion-free 
Riemannian geometry, and it can be derived by Lagrangian method. 

The torsion components can also be expressed through the stress functions by means 
of the relation 'strain-stress' and the Kroner ansatz. The direct calculation in the linear 
approximation shows us that these components are not zero, though vanishing can be 
asked for a special choice of the stress functions. In our situation, Sec. 3, vanishing occurs 
in the particular limit k/)< 1, when the modified stress functions coincide approximately 
with the Airy's and Prandtl's ones. This simply means that the torsion-less interpretation 
behind (28) is valid for (14) just at Kp <C 1. (Notice, that non-linear constitutive relation 
has to be used to re-express torsion in general situation). 

Strictly speaking, the results concerning the torsion rather imply a necessity to extend 
the framework and to pass to the Riemann-Cartan formalism, which will allow for a non- 
trivial torsion. However, the Eq.(28) itself still looks appropriate in the following sense. 
As soon as the matter is concerned with the defects sufficiently separated each from other, 
the torsion can be required as localized inside tubes confining the defects. Outside such 
tubes (i.e. for p > R cor e > 0) but not very faraway from the cores (p < R ex terior < 00), 
the Eq.(28) can approximately be valid for the purely translational defects at up <C 1. 
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Indeed, the Eq.(14) with the solutions found capture properly the required conventional 
stress functions and can be considered as a first step of an iterative scheme at k — > 0, 
provided (28) would dictate next corrections, valid outside R core . 

Recall that the stress function approach [20, 21] postulates vanishing of the Einstein 
tensor, while the torsion (i.e. the density of dislocations) is presumed 5-like to obtain the 
first order approximation solutions. This picture corresponds to teleparallel geometry. On 
the contrary, in our case the source in the Einstein equation is present, but it is smeared. 
Hence, the first order solutions appear only in the special limit k — > 0. 

It is the Hilbert-Einstein Lagrangian which provides the Einstein tensor as the L.H.S. 
of the gauge equation. However it requires a modification of the geometric picture [1, 2]. 
Practically, the use of the "coordinate" representation seems to be appropriate. Besides, 
now it is more clear, how the matter Lagrangian £g affects the linearized model, e.g., 
it leads to the limited modified stress functions instead of the increasing conventional 
ones. Although we have started with a rather special problem of modification of the 
033-component, it seems, eventually, that the ideas discussed here should be valid for the 
whole ISO(3)-mode\ also. 

As the immediate implication of the present consideration, let us mention the fol- 
lowing one. The second order corrections to the conventional screw dislocation, which 
are implemented by [1, 2], have been obtained in [7]. As far as (3) is quadratic in 
derivatives of the gauge potential, the L.H.S. of the gauge equation is of first order in 
second derivatives. Therefore, it remains the same for both perturbative steps in [7]. 
But the Lagrangian (3) is shown to be fairly inconsistent when replicating the edge 
dislocation, while the Hilbert-Einstein one removes the problem. Therefore usage of 
suggested here could influence the results of [7], because the L.H.S. of the corresponding 
Einstein-type equation is non-linear. Moreover, the use in [7] of the constitutive relation 
without the Lame constants of third order could be restrictive. 

Generally speaking, both the Lagrangians, the Hilbert-Einstein and the quadratic, can 
occur in the gauge-translational models. The matter should be to specify a proper range 
of problems for each of them. The consideration presented is to show that the Hilbert- 
Einstein Lagrangian is rather appropriate when replicating the conventional dislocations, 
but further investigations are needed to go beyond the weak field. In its turn, the La- 
grangian (3) can be more appropriate in such problems as [3-5]. Besides, the Refs.[40] 
should be mentioned where a Lagrangian quadratic in the translational field strength has 
been used in more general form to study a non-gravitational structure on space-time. 

To summarize, a non-linear Lagrangian T(3)-gauge model is proposed which admits 
at rep < 1, the conventional dislocation solutions of the linear isotropic (incompatible) 
elasticity. The calculation presented can be considered as a first order approximation. 
Incorporation of second order contributions to the gauge equation, and/or extension to 
the 75'0(3)-case to include the torsion are needed to extend the present consideration, 
and to decide about further perspectives of the 7S'(9(3)-gauging for defects modelling. 
With regard to the great attention to the T(3)-gauge models with the quadratic it 
is hopeful that the work presented would serve to a more adequate understanding of the 
problem discussed. 

A historical remark: more superficial versions of the critical content of the Sec. 2 and 
of the correction discussed in Sec. 3 can be found, respectively, in [41] and [42]. 
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